Free energy and torque for superconductors with different anisotropies of H c2 and A 
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The free energy is evaluated for a uniaxial superconductor with the anisotropy of the upper critical 
field, 7.S = H C 2, a b/Hc2,c, different from the anisotropy of the penetration depth 7a = X c /\ a i,. With 
increasing difference between jh and 7a, the equilibrium orientation of the crystal relative to the 
applied field may shift from 8 — tt/2 (8 is the angle between the field and the c axis) to lower angles 
and reach 8 = for large enough jh- These effects are expected to take place in MgB2. 
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PACS numbers: 74.60.Ec, 74.20.-z, 74.70.Ad 

It is a common practice to characterize anisotropic su- 
perconductors by a single anisotropy parameter defined 
as 7 = ^q/^c = A C /A Q (£ is the coherence length, A is 
the penetration depth, and a, c are principal directions 
of a uniaxial crystal of the interest here). The prac- 
tice emerged after the anisotropic Ginzburg-Landau (GL) 
equations were derived for arbitrary gap and Fermi sur- 
face anisotropies by Gor'kov and Melik-Barkhudarov [Q. 
Formally, this came out because in the GL domain, the 
same "mass tensor" determines the anisotropy of both £ 
(of the upper critical fields H C 2) and of A. 

At arbitrary temperatures, however, the theoretical 
approach for calculating H c2 (the position of the sec- 
ond order phase transition in high fields) has little in 
common with evaluation of A (the weak-field relation 
between the current and the vector potential), so that 
the anisotropies of these quantities are not necessarily 
the same. In fact, in materials with anisotropic Fermi 
surfaces and anisotropic gaps, not only the parameter 
7# = H C 2, a f H C 2, C may strongly depend on T, but this ra- 
tio might differ considerably from 7a = A C /A Q at low T's. 
For MgB 2 7# ~ 6 at low temperatures ||, ||, f|, || . There 
is no consensus yet on the low-T 7a. Zehetmayer et al 
use the torque data to find 7a ~ 7#; however, their pro- 
cedure implies this equality to begin with (as discussed 
in detail below) . Microscopic calculation for clean MgB2 
yields 7a » 1.1 [§, With increasing T, j H (T) de- 
creases, the calculated 7a increases till they meet at 
T = T C : 7 ff (T c )=7A(T c )«2.6. 

One of the most sensitive methods used to extract the 
anisotropy parameter 7 is to measure the torque acting 
on a superconducting crystal in the mixed state situated 
in the applied field tilted relative to the crystal axes. In 
intermediate field domain, H c \ <C H <C H C 2, the demag- 
netization shape effects are weak, and the torque density 
can be evaluated fi 01 : 
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where is the angle between the induction B and the 
crystal axis c, 



A 3 = A 2 h A c , and r\ ~ 1. This formula can be written as 
t = M x H ; since in this field domain the magnetization 
M <C H, one can disregard the distinction between B 
and H so long as the torque is concerned. It has been 
assumed in the derivation of Eq. (|l|) that the anisotropies 
of H C 2 and of the London penetration depth are the same: 
1h = 7a = 7- 

Note that Eq. (||) describes the system with the stable 
equilibrium at 6 = tt/2, i.e., the uniaxial crystal in the ex- 
ternal field positions itself so that the field is parallel to ab 
(if 7 > 1). It is also worth noting that if one applies for- 
mally the expression ([!]) with H C 2{6) — H C 2, a /e(d) having 
the anisotropy different from that of A, one may obtain 
qualitatively different angular behavior of t(0). 

Expression (0) for the torque has been derived within 
the London approach by employing the cutoff at dis- 
tances on the order of the coherence length £ where 
this approach fails; that is how the upper critical field 
Hc2 ~ 4>q /£ 2 enters the London formula. The formula, 
however, has been confirmed experimentally with a good 
accuracy as far as the angular dependencies of quantities 
involved are concerned [[To]. Uncertainties of the Lon- 
don approach are incorporated in the parameter r\ ~ 1; 
discussion of those can be found, e.g., in Ref. [ll]. 

Below, the free energy of the mixed state and the 
torque are derived for the general case of 7// 7^ 7a- The 
torque expression is shown to acquire new terms which 
describe a more complicated behavior as compared to 
that of Eq. (Q). This implies that the analysis of the 
torque data might be misleading if Eq. (Q) is employed 
to materials like MgB 2 §, [TJ|, [jjj. 

Let us start with the London expression for the free 
energy valid for intermediate fields H c i <C H <C H C 2 
along z tilted with respect to the c crystal axis over the 
angle toward the crystal direction a pi : 
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7a > m 2 a m c 



e(0) = J sin + 7 2 cos 2 , 



(2) 



here, m zz — m a sin 2 0+m c cos 2 0, m c /m a 
1 for uniaxial crystals, and G form the reciprocal vortex 
lattice. The summation is extended over all nonzero G. 
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As usual we evaluate the sum here by replacing it 
with an integral over the reciprocal plane: — > 
(4v/i^B)fdG,dG v : 
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G y , and we use polar 



where g x = y 7 ? 
coordinates: g x — g sirup, g y = g cos tp. 

To determine the limits of integration over g, one can 
recall that in fields B <C H C 2 the vortex lattice struc- 
ture is fixed by the anisotropy parameter 7a 14 . In 



tintermediate fields, for all possible equilibrium London 
structures, the distance in the reciprocal space g from 
the origin to the nearest neighbors is given by pa] 
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Therefore, one can take go (multiplied by a number of 
order unity) as the lower limit in the logarithmically di- 
vergent integral over g. 

The upper limit in this integral is affected by the form 
of the vortex core. To determine the core shape we note 
the microscopic evaluation of H C 2 (see Ref. |^) shows that 
with a good accuracy the angular dependence of H C 2 is 
given by the standard GL form at any T: 
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This can be written as 
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where £ XjV are semiaxes of the elliptical core. We stress 

that the "masses" are different from rriik which 

determine the anisotropy of A; in particular, fj, a = 
-2/3 4/3 , -2/3 4/3 

Tff 3 Mc = 1h , whereas m a = 7 A m c = 7/ . 

Equation (Q) gives maximum values of G^ and G y : 
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Thus, the domain of integration in the G plane is bound 
by an ellipse 
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In other words, at a given ip, the upper limit in the inte- 
gral over g is 



9m{<p) = 



£\//3 2 cos 2 tp + sin 2 ip 
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and we obtain: 
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Note that /? = 1 for coinciding anisotropics of H C 2 
and A; the integral J(j3) over 95 then vanishes and we 
have the standard expression for the energy. To evaluate 
</(/?), we observe that the integral for dJ/d(3 is a ratio- 
nal function of cos 2 <p> and therefore can be calculated 
by going to the complex plane with the help of residues: 
dJ/d/3 = 4tt/(/3 + 1). Since J(l) = we obtain: 



Thus, we have 
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To write explicitly the angular dependence of F, it is 
convenient to use the angular functions 



Ox.h(0) =ex, H {P)hx, H 



(14) 



where £\,h{@) is defined in Eq. (^) with corresponding 
7's. In terms of these functions, f3 = Oh/Q\ and 
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The torque density follows 
cj) B 
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where e 
because 



2.718. 



The torque is zero at 
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for both 9a and Oh- In the standard case of jh = 
7a = 7, Eq. ( p^[ ) reduces to the result (|l|) if we set 
77 = 7rv3/e ~ 2 . We then can rewrite the torque in the 
form: 
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FIG. 1: The torque in units of cj>oB /32ty 2 X^ b versus angle 
< 6 < tt/2 for 7 A = 2.2, -y H = 3, and 4e 2 r)H c2 , c /B = 100. 



Since Q' x H < 0, the second contribution to the torque 
( |l6| ) is negative whereas the first one is positive. The 
positive torque implies that the system energy decreases 
with increasing 9, as in the case of jh = "fx for which 
9 = tt/2 is the stable equilibrium. 

The competing roles of this two contributions can 
be demonstrated by considering stability of equilibrium 
states at 9 = and 6 — tt/2. To do this one notes that 
6(0) = 1, 6(tt/2) = l/ 7 , 9'(0) = 9'(7r/2) = 0, and 



e (o) = - 
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(for both 0\ and 0#). Then, one obtains: 
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where the constant positive prefactor is omitted since we 
are interested only in the sign of F . Clearly, 9 = tt/2 
corresponds to the stable equilibrium for — 7a- In 
the general case, however, there is no such a clear-cut 
result: for a fixed 7a and large enough 7#, 9 = tt/2 
may become unstable. E.g., for 7a = 1 and jjj > 1, 
f"(tt/2) < whereas F" (0) > . 

To illustrate how the angular dependence of the torque 
varies with anisotropies of H C 2 and A, we evaluate nu- 
merically the torque density (fLq) for parameters in the 
range of those for MgB2. Fig. [jjshows t(9) for 7a = 2.2 
and 7# = 3, the values expected for temperatures some- 
what below T c . Qualitatively, the dependence is stan- 
dard; the torque is positive in the whole angular domain, 
i.e., 9 = tt/2 is the stable equilibrium. 

With decreasing T, of MgB2 increases whereas 7a 
decreases. In Fig. |^ the torque ( |l8| ) is plotted for 7a = 2, 
jH = 5 (the upper curve) and for 7a = 1.7, 7# = 5.3 (the 
lower curve). These values correspond roughly to 0.7 and 
0.6 T c according to Ref. 0. Clearly, 9 = tt/2 as well as 
9 = are unstable; the stable equilibrium is shifted to 
< 9 < tt/2. Interestingly enough, Angst et al. observe 
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FIG. 2: The same as in Fig.jl]. The upper curve is calculated 
with Eq. (|l^) for 7a = 2 and jh = 5; at the lower curve 
7a = 1.7 and = 5.3 . 




FIG. 3: The same as in Fig.jl], but 7a = 1.1 and jh = 6. 



that a strong "peak-effect-like" irreversibility develops in 
the torque data at » 77° at T = 15 K and H = 7.5 T [0. 
In layered materials this effect is commonly interpreted 
as manifestation of the "intrinsic pinning" in the small 
vicinity of the equilibrium orientation at 9 = tt/2. From 
the point of view of this article, the peak should move to 
a position of the stable equilibrium, i.e., to lower angles. 
Moreover, the data show a negative torque above this 
angle (77° < 6 < 90°) in agreement with Fig. [| Still, 
interpretation of these data within the London model 
should not be taken too seriously: the applied field H = 
7.5 T exceeds H c2yC (0) and in the ^-domain where H < 
Hc2{0): the London model should not be trusted. 

Finally, we plot in Fig. || the torque density for pa- 
rameters which correspond to low temperatures, where 
7a ~ 1.1 and 7# sa 6. The torque is negative for all 
angles implying the stable equilibrium at 9 = 0. 

Physically, the large low temperature anisotropy of H C 2 
in MgB2 is caused by the large superconducting gap on 
the nearly two-dimensional sheets of the Fermi surface 
of this material 0, || ||. With increasing T, ther- 
mal mixing with the states at the three-dimensional part 
of the Fermi surface suppresses this anisotropy to about 
1h{T c ) fa 2.6. The anisotropy of the London A (or of the 
superfluid density) at T = of clean materials does not 
depend on the gap at all ( "Galilean invariance of the su- 
perflow" ) and therefore is determined by the whole Fermi 
surface, i.e., it is weak for MgB2 (see discussions in Ref. 
g, 0). The calculation § shows that 7A (T = 0) w 1.1 
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and grows to w 2.6 as T — > T c . 

Thus, different gaps at different Fermi surface pieces 
or keneraffy, anisotropic gaps on anisotropic Fermi sur- 
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( 

face) may lead to profound macroscopic consequences 
such as those considered above. Certainly, the equilib- 
rium magnetization and the flux lattice structure should 
be strongly affected by differences in anisotropics of the 
upper critical field and of the London penetration depth. 
Occurrence of a peak-effect near the field orientation 
other than H \\ ab is an example of peculiar dynamic 
phenomena which call for further study. 
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